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1 Introduction 

More than a decade ago, Grover [5j discovered a quantum algorithm for 
searching unsorted list runs in 0{^fN). The algorithm is important as the 
best that can be attained classically with an unsorted list is in the order 
of 0{N). Since the publication of the work, there have been various im- 
provements made to Grover 's algorithm. Some work ( [2] , [3] , [9] , [TO] ) extends 
the result to other search-based algorithms, while others ([3], [7], [8], [13]) fo- 
cusses on the nature of the unitary transformation involved in the search 
algorithm. The latter research direction is to be expected as in the words of 
Grover ^ , the core of the algorithm is the " design of the unitary evolution 
of the system" . 

Though Grover's algorithm clearly outperforms its classical counterpart, 
it is still not 'fast' enough when applied to NP-complete problems. A 
straightforward application to the Travelling Salesman Problem, for exam- 
ple, following the approach in [3], results in an 0(\/nI) performance, where 
n is the number of cities. The result is still not computationally tractable 
for large values of n. This is surprising, as due to the inherent parallelism 
in quantum computation, one would expect a better result. Hence, in this 
paper, instead on focussing on the design of the unitary search evolution, we 
focus instead on the possibility of embedding a quantum search algorithm 
within a classical binary search framework. The result appears promising: 
taking full advantage of quantum parallelism, we show that it may actually 
be possible to search an unstructured list in 0{lg{N)), provided we are will- 
ing to restart the quantum search multiple times with a different sequence 
of qubits and perform a series of measurements at the end of each. 

The gist of the idea in our research approach comes from a classical 
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algorithm within Computer Graphics and the Computational Geometry: 
ray tracing. In ray tracing, a ray is casted into a scene comprising of a 
number of objects. A brute force approach would test the ray against each 
and every object in the scene, a process that scales linearly with the number 
and the complexity of the objects. A smarter approach requires that the 
scene be partitioned into regions, and the bounding volume of each region 
be approximated. For each region, the ray would then be tested first against 
the bounding volume approximation. If the test is positive, then only would 
the ray be tested against each object within the region. Can a similar scheme 
be invented for Grover's algorithm? We would like to partition the list of 
items being searched into sublists. For each sublist, before we even perform 
repetitive amplification, we first check as to whether or not the desired item 
is somewhere within it. If the item is not there, we move on to the next 
sublist, and repeat the same query. We elaborate on this idea in the next 
section. 



2 Binary Subdivision Algorithm 

Our algorithm can be intuitively understood in classical terms as follows: 
We wish to search for an item, t, in an unsorted list, L. Suppose we are 
given a 'magic' function, F, that given L as input, instantaneously provide 
an output that indicates whether or not t is within the list. F does not 
pinpoint the location of t; it merely says whether or not t is somewhere 
within the input list. We can then search using the following procedure: 
Search (L) 

1. Divide L into two sublists: L\ and L2. 

2. Using F, check whether Li contains t. If yes, call Search(Li). If not, 
call Search(L-2). 

Assuming that t exists in L, the recursive procedure will lead to t. We 
claim that F can be implemented as a quantum function, and that the 
overall procedure can be run classically. 

As in Grover's algorithm, we start with 2 registers: n qubits in the 
first, and 1 qubit in the second. Note that the n qubits in the first register 
represents N = 2^ numbers. In Grover's algorithm, the first register is 
initialized to be in the state \Q) = |0)®", and the Hamadard operator, i/®", 
applied to it. The result is a linear combination of all 2" computational 
basis states \S) = -^X^^q^K)- The second register is initialized to 
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and after a Hamadard transformation, comes to be in the state |— ). Its sole 
purpose is to serve as the second qubit for the It oracle operator. Repeated 
application of Grover operator, G = —HIoHIt, is performed to increase the 
probability of obtaining the desired state, \t), when a measurement is made 
upon the first register. The If in the operator inverts the phase of \t). The 
Iq operator inverts |0). Altogether, Grover operator works by amplifying 
the phase of the desired item, 

Instead of performing a quantum search on the entire list of items within 
\S), we wish to search only where \t) might be. Our approach requires two 
main features: 

1. the ability to determine whether or not a state \S) contains \t). The 
function F described in the first paragraph of this section has this 
ability. 

2. the ability to partition the list of items being searched into sublists. 

The first feature is implemented by first running the oracle It in such a 

way that its second register is |1) for a target item, and |0) otherwise. As- 
suming that at most there can be only one target item, the resulting state 
of the second register will then be either |0) or yl^^l/iV |0) + 1/y/N 
The distance between the 2 possible results decreases exponentially with 
increasing value of N . We show however how a series of non-unitary mea- 
surement operators with non-negligible probabilities of success can be used 
to differentiate between the two possible states. 

As for the second feature, we note that the items in Grover's search 
do not actually exist as in a classical, conventional list. Instead the list is 
actually a single state, \S), that encodes a superposition of possible item 
values. We can however perform a partitioning by fixing qubits within \Q) 
and appropriately modifying the form of the Hamadard transform applied 
to it. As a specific example, consider the leftmost qubit, that is the one with 
the most significant digit. Half of the items in the state \S) that result after 
a H'^'^ operation on \Q) will start with |0) and the other half with But 
if we fix the leftmost qubit in \Q) to be |0) (or |1)) and perform a I®II®'^~^ 
instead, then the state IS') will only consists of items that start with |0) (or 

A more elaborate, formal elaboration of the proposed algorithm is then 
as follows: Let the initial state of the first register be \Q) = jgo) \<li) ■■■ \<ln-i)- 
Assume that we are now determining the value of the k*^ qubit within the 
first register, where Q < k < n. Let h he a, binary variable that contains 
the state value (|0) or to be tried out for qj^. Let \Qq) be the sublist of 
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\Q) for which the qubit values have already been determined. The general 
form of \Q) is then as follows: {\Qo) \qk) kfc+i) ■•■ \Qn-i))- Define as usual 
the function, / : |0) , |iV — 1) ^ 0, 1, that recognizes the solution: 

I otherwise 
Our algorithm then proceeds as follows: 

1. Initialize \Qo) to be a null list (|)), and set A; = 0, and 6 = 0. Initialize 
as well the second register to |0). 

2. Repeat the following steps: 

(a) Prepare the initial state \Q) = \Qo) \b) |0) |0) ... |0). 

(b) Perform I^'' (g) H^i^''''^ on \Q): 



im ^ H^in-k)^^Q^^ 1^,^ |o) ... |0)) = IQo) |6) = IV^) 

(2) 

(c) Perform It on l'^) |0): 

m |o) = (^Qo) \b) ^g=^ 

_ 1^ . S-:o'"K) |oe/(|Qo)|6)K))) 

(d) Consider the 2 possible results when the second register is mea- 
sured: 

• If the target item is not within {ip), then the result of the 
observation will be \y) = |0). Otherwise, it will be 



where N = 2^. Assume for now that it is possible to dif- 
ferentiate between |x) and \y). Hence, we can proceed as 
follows: 

— Suppose the target item is not within \tp). Then if b == 
0, we set = 1 and repeat step 2. 
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— If the target item is within 1-0), we then append \qk) to 
IQo) - its state value considered fixed -, and proceed on 
by repeating step 2 with 6 = and k = k + 1. 

The complexity of the above algorithm, in terms of the number of queries 
to the oracle, assuming that we perform recursive subdivision without actu- 
ally doing any amplitude amplification, is obviously lg{2^). 



3 Differentiating between \x) and \y) 

We need to diff'erentiate between 



- 1 1 

|0> + ^|1) (6) 



and 

\y) = |o) (7) 

The dot product between the two however approaches 1 with increasing 
N , that is ■\/ {N — 1)/N — t- 1 when N — )• oo. Nevertheless, if it is possible 
to apply a transformation using the following matrix, we would be able to 
tell the two apart: 

Note that: 

and 



. VN ^ 

D\x) = \l) (9) 
D\y) = \0) (10) 



Unfortunately, D is not a unitary matrix. It is however possible to apply 
a non-unitary operator as a quantum measurement operator, as proposed in 
the work by Terashima and Ueda [12] . We need to convert D into an operator 
that is implementable as a quantum operator. To do this we perform the 
following transformation steps: 

1. We first factor D using Singular Value Decomposition (SVD): 

D = QVR^ (11) 

Q and R are unitary matrices, while y is a diagonal matrix. Q and R 
can be implemented easily within a unitary quantum framework. We 
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are left with V, a2x2 diagonal matrix that comes with one very large 
and one very small diagonal value. For example, if n = 10, V would 
be: 

/ 45.249308037472204 ^ ^21 



0.707193134832027 

And if n = 30, V would be: 

/ 46340.95000644678 , , . 

~ \ 0.70710678127 ' ^ ' 

2. We now factor V into its roots, such that the first diagonal element 
becomes a value less than 2. The power of the root depends on the 
size of the first diagonal element. If we take the 16*^^ root of V, for 
example, we obtain 

V = 14/16^^1/16-^1/16 (14) 

where Vi/ig is applied 16 times. 

In the case of n = 30, Vi/iq would be as follows: 

_ / 1.957144124161160 \ , . 

^/^^ ~ \ 0.978572062094820 J ^ ' 

3. We then normalize the diagonal elements of 1^/16 • L^t Mq be the 
resulting matrix. In the case of n = 30, Vi/i6 matrix would then 
become 

_ / 0.894427190997313 ^ ri6l 

° ~ V 0.447213595505164 ) *■ ^ 

The sequence of transformation to be applied to either |a;) or \y) is then 
in the following order: 

1. i?t 



2. M0...M0 {v times where v indicates the power of the root of V applied 
to 

3. Q 



to obtain Vi/^) 
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We prepare the matrices used in the above sequence for each of the 
possible values of n encountered in the quantum binary search in Section 2. 

Note that the measurement operator complementary to Mq is computed 
as follows [12j: 



Ml = ^1 - Mq^Mo (17) 
The effect of Mq on a state is computed as follows: 

(18) 



^ (V'lMjMol^) 
3.1 Example 

We consider as a complete numerical example, the case when n = 20. Sup- 
pose the target item is within the state in the first register. The state of the 
second register will then be: 

X = 0.999999523162728 |0) + 0.000976562500000 |1) (19) 

The matrix D that we want is: 

/I -1023.999511718634 \ 
^ ~ V 1024.0 J ^ 

Performing SVD, we obtain D = QVR^ where 

/ -0.707106781186547 -0.707106781186547 \ 
^ ~ y 0.707106781186547 -0.707106781186547 J ^ ^ 

f 1448.154515236507 \ 

~ I 0.707106865480 j ^ ' 



and 



Now, 



_ , -0.000488281308208 -0.999999880790675 . 
~ ^ 0.999999880790675 -0.000488281308208 ' ^ ' 



1/16 _ / 1.575980833365910 \ 

» 0.978572069378633 ' ^ ' 



Normalizing the matrix (along the diagonal), we obtain Mq, a non- 
unitary measurement operator: 

_ / 0.849549077650853 \ 

° ~ i 0.527509587270776 j ^ ' 
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With the matrices above prepared (in an offline process), we follow the 
proposed sequence of steps. We start off by applying the unitary operator 
R to obtain: 

\d) = 0.000488281308208 |0) - 0.999999880790675 |1) (26) 

We go now into the non-unitary phase, where Mq is to be applied 16 
times. Consider the effect of just the first iteration. Note that the success 
probability in applying Mq is 0.278266470393446. After applying Mq, state 
becomes: 

\d) = 0.000786372165026 |0) - 0.999999690809361 11) (27) 

We assume for now (in this paper) that Mq is successful throughout the 
non-unitary phase. Note that the success probability in applying Mq in this 
particular example is in the range from about 0.28 to about 0.40. 

At the end of the 16 applications of Mq, the state becomes: 

\d) = 0.707106781186547 |0) - 0.707106781186548 |1) (28) 
Applying Q, the state finally becomes: 

Id) =0|0) + 1|1) (29) 

4 Notes 

We make a few notes with regard to the method proposed in this paper. 

1. Quantum search using binary trees has been considered in other works, 
for example [9\- However, the binary tree search proposed in this paper 
does not assume sorted data. 

2. An alternative approach for differentiating between \x) and \y) is de- 
scribed by Ohya and Volovich in [11]. They propose the use of what 
they call a chaotic amplifier to separate the two states. The method 
proposed in this paper, however, should be easier to understand and 
implement. 

3. While it has been proven before that a black-box quantum query al- 
gorithm cannot solve NP problems in o(2"/^) [1], the approach in this 
paper incorporates non-unitary model, and hence it does not violate 
the result in [1]. 
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4. The main catch in the proposed method is that in the state differen- 
tiation module, we have not considered the case when Mi instead of 
Mq was apphed during the non-unitary phase, resulting in the inter- 
mediate result \d') rather than \d). In some cases, \d') is far enough 
from \d) so as to be easily distinguishable. In general, some sort of 
statistical scheme would be needed to ensure that the final result used 
is one obtained without or with minimal application of Mi . 
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